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Meridional flow profile measurements with SOHO/MDI 
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We present meridional flow measurements of the Sun using a novel helioseismic approach for analyzing SOHO/MDI data 
in order to push the current limits in radial depth. Analyzing three consecutive months of data during solar minimum, 
we find that the meridional flow is as expected poleward in the upper convection zone, turns equatorward at a depth of 
around 40 Mm 0.95 Rq), and possibly changes direction again in the lower convection zone. This may indicate two 
meridional circulation cells in each hemisphere, one beneath the other. 
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1 Introduction 

The first observations of the Sun's meridiona l flow were 
made independ ently by Duvall (Il979h . Howard ( 19791) . and 



Beckers (Il979l) . who found a poleward drift of r^lOm/s 



at the solar surface. By mass conservation, this meridional 
flow has to change direction somewhere in the solar interior, 
which may happen in the lower regions of the convection 
zone. 

The meridional circulation plays an important role in 
many solar models. In mean-field models of angular mo- 
mentum transport, the meridional circulation is responsi- 
ble for how differen t ial rotation is established and main - 
tained (Rudiger'l989': lThompson et al.l2003l:lMieschl2005h . 
Anisotropic Reynolds stresses, arising from random turbu- 
lent fluctuations of the plasma motions, can induce a merid- 
ional circulation, which in turn drives the differential rota- 
tion. In numerical simulations, Brun & Toorm'e ( 2002.) find 
that the differential rotation profile is established by turbu- 
lent convection, and that the Reynolds stresses play a crucial 
role in transporting meridional momentum toward the equa- 
tor. The meridional circulatio n is also an importan t ingredi- 
ent in some dynamo models (ICharb onneau 2005'). In flux- 



transport models of the solar dynamo (Babcock-Leighton 
models) the meridional flow is needed to transport the sur- 
face poloidal magnetic field to the bottom of the convection 
zone, where it is then converted into a toroi dal field by rota- 
tional shear ( Dikpati & Charbonneaull9"99 ). 

Using helioseismic techniques, the meridional flow is 
usually observed in the upper layers of the solar convec- 
tion zone. Quite good spatial resolution is achieved in this 



needs to be large. No return flow has conclusively be en ob- 
served to date. Using time-distance helioseismology, [Giles 
(Il999h measured the radial profile of the meridional circu- 
lation down to the bottom of the convection zone and, with 
the added constraint that the meridional circulation has to 
be contained within the convection zone, inferred a return 
flow at around 0.8 solar radii. Applying Hankel analysis. 



Braun & FanI (119981) saw a tentative return flow at around 



40 Mm depth (bu t also c onsistent with zero within the error 
margin), whereas Duvall & KosovichevI ( 1200 lb observed no 
return flow. 



2 Data and analysis 

We analyz e velocity data from the Michelson-Doppler Im- 
ager (MDI.I Scherrer et al.lll995h on board flie Solar an d He- 
liospheric Observatory (SOHO. iDomingo et alJll995b . The 
observations were taken during three consecutive Carring- 
ton Rotations: CR 1 922 ( 1 997 April 24 - May 2 1 ), CR 1 923 
(1997 May 21 - June 17), and CR 1924 (1997 June 17 - 
July 13). During this time, the Sun was at the minimum 
of its activity cycle. We use tracked data along the Center 
MeridiarQ, obtained for the SOI Dynamics Campaign us- 
ing high-resolution, full-disk Dopplergrams. The cadence is 
60s. 



2.1 Dispersion relation along the center meridian 

As our aim is to measure the meridional flow, which is in 
North-South direction, we restrict ourselves to analyzing 



regime, in particular with ring-diagram analysis (e. g.. lGonzaIezTO tMte^ftl<^.| i" the same direction. The data are co- 
1999h . Deeper down in the convection zone the horizontal added along latitude, resultin g in a 2D array with one time 
resolution is poorer and to probe there the viewing angle and one spatial coordinate (see|Mitra-Kraev et al.|200^. The| 
spatial coordinate 9 (latitude) is equally spaced in angle 
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Fig. 1 The difference between the negative and the positive power spectrum, Pneg — Ppos- The frequency shift is visible 
along the f- and p-mode ridges. Dark dots are for Pnog > ^pos, while lighter dots are for Pncg < Ppos- Northward flow is 
present if the dark dots are shifted to higher frequencies and lighter dots to lower ones, while for Southward flow lighter 
dots are at higher frequencies than dark dots. The size of the dots is proportional to the absolute power The black lines 
indicate the lower ray-turning points. 



(spacing ~ 0.125°). As the meridional flow is oppo- 
site in each hemisphere and vanishes around the equator, 
flie range +20° to +60° is used for the Northern, and -60° 
to —20° for the Southern hemisphere. For each hemisphere 
we obtain a power spectrum 



where 

f{>^, ke) 



-i{2-iTut+ke9) 



is the Fourier transform of the data /(t, 6). The power spec- 
trum has symmetries such that P{\v\ > 0, |^| > 0) = 
< 0, \l\ < 0) = Ppos, which we define as the 'pos- 
itive', and P{\iy\ > 0, |/| < 0) = P{\iy\ < 0, |/| > 0) = 
^ncg, which we define as the 'negative' power spectrum. 
The difference between the positive and the negative power 
spectrum can be attributed to flows. 

The approach of a ID Fourier transform along the center 
meridian is comparable to a spherical harmonic treatment 
with only con sidering rn = Q modes, which was used by 
Krieger et al.1 tOOl\) for meridional flow measurements. As 
kg corresponds to the spherical harmonic degree /, we use 
the notation kg = I from now on. With the observing length 
Tobs ~ 27 days and the latitudinal range Grange — 40°, 



the resolutions are given by dv = l/Pobs ~ 0.4 /iHz and 

dl = 27r/0,a„gc « 9. 

2.2 The effect of flows on the p-mode ridges 

For slow flows, which allow a linear approximation, the fre- 
quency splitting is given by Auj ^ kh ■ U' (an expression 
which is also used in ring-diagram analysis) where kh is the 
horizontal wave number and the velocity U' is the radially- 
dependent horizontal flow U (r) convolved with a function 
describing the medium the pressure wave passes through. 
For measured Aui as a function of kh, U (r) can be retrieved 
by an inversion. 

In our case kh, U' and U (r) are ID variables, describing 
the flow perpendicular to the solar radius in the North-South 
direction. With Aw = 27rAi/ and kh — 1/Rq we obtain for 
the frequency shift in velocity units 



(1) 



where Rq = 6.96 x 10^ m is the solar radius. The fre- 
quency shift IS.V is the shift by which, for a given I value, a 
p-mode ridge is displaced from its unperturbed position due 
to a flow. A Northward flow moves p-mode ridges along 
/, such that the frequency of p-mode ridges in the positive 
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Fig. 2 The observed frequency shifts in the Northern (top) and Southern (bottom) hemispheres. 



power spectrum is reduced by Ai^, while the frequency of 
p-mode ridges in the negative power spectrum is enhanced 
by At/. 

Figure [T] shows the difference between the negative and 



the positive power spectrum, Pn, 



Ppos, for the South- 



ern hemisphere of CR 1922-1924. The size of the dots is 
Hnear with absolute spectral power. For better readability, 
all power smaller than a lower cutoff, which was chosen 
to decay exponentially along I and has a Lorentzian shape 
as a function of v around 3.3 mHz, was set to zero. The f- 
and p-mode ridges are clearly visible due to the frequency 
shift. The f mode as well as pl-p8 modes are labeled. Dark 
dots are for P, 



neg 



^pos > 0, while lighter dots repre- 
sent Pneg — ^pos < 0. This means a lighter dot being at 
a slightly larger frequency than its darker counterpart in- 
dicates a Southward flow, whereas a dark dot sitting at a 
larger frequency indicates a Northward flow. The straight 
black Unes indicate the lower ray-turning points. 

2.3 Measuring tlie frequency sliift 

In order to obtain the frequency shift Ai^, we cross-correlate 
Pnog and Ppos for each I over a z^-range covering each p- 
mode ridge. The criterion for accepting a measurement is a 
correlation coefficient of 0.35 or more. 

Figure |2] shows the measured frequency shift U'{v/l) 
for the Northern (top panel) and Southern (bottom panel) 
hemispheres obtained for each power spectrum of the three 



Carrington rotations. The spread (—200 to +400 m/s) of 
the individual frequency shifts, represented by crosses and 
color coded with the corresponding p-mode, exceeds the 
plotting range. The middle black curve was obtained by 
smoothing the individual measurements twice with a bin 
size of 101 points. For the first smoothing, the error is given 
by the standard deviation of the points. The error margins 
in Figure 12] (upper and lower lines) are then given by error 
propagation in the second smoothing, assuming that the ini- 
tial errors given by the standard deviation of the points are 
independent. However, as the initial errors are clearly not 
independent, this treatment undoubtedly underestimates the 
errors somewhat. 



2.4 Asymptotic inversion 

We invert the smoothed curve U' using an asymptotic in- 
version (see, e.g.. Christensen-Dalsgaard et al..,1990.) . The 
inversion is given by 



U{t) 



^2a{r) d r^'') 



1 



dlnr Ja(RQ) V a2(r) - 



U'{w) dw 



with w = 2'Kv/l, r the solar radius, a{r) = c{r)/r, and c(r) 

the sou nd-speed profile given by Model S of Christensen-Dalsgaard et aL^ 
1996h . The integral ranges from the solar surface to the 



lower turning point of the ray path. In order to solve the 
integral, which has an integrable singularity at the upper 
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Fig. 3 The meridional flow profile for the Northern and Southern hemisphere obtained by asymptotic inversion. 



boundary, the numerical scheme described in Appendix II 
ofiChrist ensen-D als gaard et al, ( 1989) was applied. 

Figure [5] displays the obtained meridional flow profile 
from the asymptotic inversion. Displayed is the Northward 
flow for both the Northern and Southern hemispheres. The 
error region (shaded area) was obtained by Monte-Carlo 
simulations. The inversion profile was twice smoothed. 



3 Discussion and conclusions 

We analyzed high-resolution full-disk velocity data from 
SOHO/MDI along the center meridian covering three Car- 
rington rotations from April to June 1997. Using a Fourier 
transform and measuring the frequency shift of the p-mode 
ridges in the power spectrum, we are able to obtain a merid- 
ional flow profile as a function of radius/depth. Figure [3] 
shows that in the Northern as well as the Southern hemi- 
sphere the flow is poleward close to the solar surface and 
turns equatorward at around 40 M m. Such a re turn flow was 
previously tentatively observed by lBraun & Fa n (1998), an- 
alyzing data from the same period (but only for 8 days in- 
stead of 3 months). The flow may again turn poleward in 
deeper regions indicating a counter cell, although the data is 
also consistent with zero flow at lower regions. Such co unter 



of the inclination of the solar rotation axis (IBeck & Giles 
2005h . 



cells are consistent with recent numerical simulations tMiesch 



2007h . The near-surface flow being somewhat larger in the 
Southern than in the Northern hemisphere may be because 
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